Search for: All records

ABSTRACT Fulfilling the rich promise of rapid advances in time-domain astronomy is only possible through confronting our observations with physical models and extracting the parameters that best describe what we see. Here, we introduce redback; a Bayesian inference software package for electromagnetic transients. redback provides an object-orientated python interface to over 12 different samplers and over 100 different models for kilonovae, supernovae, gamma-ray burst afterglows, tidal disruption events, engine-driven transients among other explosive transients. The models range in complexity from simple analytical and semi-analytical models to surrogates built upon numerical simulations accelerated via machine learning. redback also provides a simple interface for downloading and processing data from various catalogues such as Swift and FINK. The software can also serve as an engine to simulate transients for telescopes such as the Zwicky Transient Facility and Vera Rubin with realistic cadences, limiting magnitudes, and sky coverage or a hypothetical user-constructed survey or a generic transient for target-of-opportunity observations with different telescopes. As a demonstration of its capabilities, we show how redback can be used to jointly fit the spectrum and photometry of a kilonova, enabling a more powerful, holistic probe into the properties of a transient. We also showcase general examples of how redback can be used as a tool to simulate transients for realistic surveys, fit models to real, simulated, or private data, multimessenger inference with gravitational waves, and serve as an end-to-end software toolkit for parameter estimation and interpreting the nature of electromagnetic transients. 

ABSTRACT We identify the progenitor star of SN 2023ixf in Messier 101 using Keck/NIRC2 adaptive optics imaging and pre-explosion Hubble Space Telescope (HST)/Advanced Camera for Surveys (ACS) images. The supernova, localized with diffraction spikes and high-precision astrometry, unambiguously coincides with a progenitor candidate of $$m_\text{F814W}=24.87\pm 0.05$$ (AB). Given its reported infrared excess and semiregular variability, we fit a time-dependent spectral energy distribution (SED) model of a dusty red supergiant (RSG) to a combined data set of HST optical, ground-based near-infrared, and Spitzer Infrared Array Camera (IRAC) [3.6], [4.5] photometry. The progenitor resembles an RSG of $$T_\text{eff}=3488\pm 39$$ K and $$\log (L/\mathrm{L}_\odot)=5.15\pm 0.02$$, with a $$0.13\pm 0.01$$ dex ($$31.1\pm 1.7$$ per cent) luminosity variation at a period of $$P=1144.7\pm 4.8$$ d, obscured by a dusty envelope of $$\tau =2.92\pm 0.02$$ at $$1\, \mu \text{m}$$ in optical depth (or $$A_\text{V}=8.43\pm 0.11$$ mag). The signatures match a post-main-sequence star of $$18.2_{-0.6}^{+1.3}\, \mathrm{M}_\odot$$ in zero-age main-sequence mass, among the most massive SN II progenitor, with a pulsation-enhanced mass-loss rate of $$\dot{M}=(4.32\pm 0.26)\times 10^{-4} \, \mathrm{M}_\odot \, \text{yr}^{-1}$$. The dense and confined circumstellar material is ejected during the last episode of radial pulsation before the explosion. Notably, we find strong evidence for variations of $$\tau$$ or $$T_\text{eff}$$ along with luminosity, a necessary assumption to reproduce the wavelength-dependent variability, which implies periodic dust sublimation and condensation. Given the observed SED, partial dust obscuration remains possible, but any unobstructed binary companion over $$5.6\, \mathrm{ M}_\odot$$ can be ruled out. 

ABSTRACT Hydrogen-poor superluminous supernovae (SLSNe) are among the most energetic explosions in the universe, reaching luminosities up to 100 times greater than those of normal supernovae. This paper presents the largest compilation of SLSN photospheric spectra to date, encompassing data from the advanced Public ESO Spectroscopic Survey of Transient Objects (ePESSTO+), the Finding Luminous and Exotic Extragalactic Transients (FLEET) search, and all published spectra up to December 2022. The data set includes a total of 974 spectra of 234 SLSNe. By constructing average phase binned spectra, we find SLSNe initially exhibit high temperatures (10 000–11 000 K), with blue continua and weak lines. A rapid transformation follows, as temperatures drop to 5000–6000 K by 40 d post-peak, leading to stronger P-Cygni features. Variance within the data set is slightly reduced when defining the phase of spectra relative to explosion, rather than peak, and normalising to the population’s median e-folding decline time. Principal Component Analysis (PCA) supports this, requiring fewer components to explain the same level of variation when binning data by scaled days from explosion, suggesting a more homogeneous grouping. Using PCA and K-means clustering, we identify outlying objects with unusual spectroscopic evolution and evidence for energy input from interaction, but find no support for groupings of two or more statistically significant subpopulations. We find Fe ii  $$\lambda$$5169 line velocities closely track the radius implied from blackbody fits, indicating formation near the photosphere. We also confirm a correlation between velocity and velocity gradient, which can be explained if all SLSNe are in homologous expansion but with different scale velocities. This behaviour aligns with expectations for an internal powering mechanism. 

Abstract We present a sample of 34 normal Type II supernovae (SNe II) detected with the Zwicky Transient Facility, with multiband UV light curves starting att≤ 4 days after explosion, and X-ray observations. We characterize the early UV-optical color, provide empirical host-extinction corrections, and show that thet> 2 day UV-optical colors and the blackbody evolution of the sample are consistent with shock cooling (SC) regardless of the presence of “flash ionization” features. We present a framework for fitting SC models that can reproduce the parameters of a set of multigroup simulations up to 20% in radius and velocity. Observations of 15 SNe II are well fit by models with breakout radii <10¹⁴cm. Eighteen SNe are typically more luminous, with observations att≥ 1 day that are better fit by a model with a large >10¹⁴cm breakout radius. However, these fits predict an early rise during the first day that is too slow. We suggest that these large-breakout events are explosions of stars with an inflated envelope or with confined circumstellar material (CSM). Using the X-ray data, we derive constraints on the extended (∼10¹⁵cm) CSM density independent of spectral modeling and find that most SN II progenitors lose $\dot{M}<{10}^{-4}{M}_{\odot }{\mathrm{yr}}^{-1}$up to a few years before explosion. We show that the overall observed breakout radius distribution is skewed to higher radii due to a luminosity bias. We argue that the ${66}_{-22}^{+11}\%$of red supergiants (RSGs) explode as SNe II with breakout radii consistent with the observed distribution of RSGs, with a tail extending to large radii, likely due to the presence of CSM. 

Abstract We present panchromatic optical + near-infrared (NIR) + mid-infrared (MIR) observations of the intermediate-luminosity Type Iax supernova (SN Iax) 2024pxl and the extremely low-luminosity SN Iax 2024vjm. JWST observations provide unprecedented MIR spectroscopy of SN Iax, spanning from +11 to +42 day past maximum light. We detect forbidden emission lines in the MIR at these early times while the optical and NIR are dominated by permitted lines with an absorption component. Panchromatic spectra at early times can thus simultaneously show nebular and photospheric lines, probing both inner and outer layers of the ejecta. We identify spectral lines not seen before in SN Iax, including [Mgii] 4.76μm, [Mgii] 9.71μm, [Neii] 12.81μm, and isolated Oi2.76μm that traces unburned material. Forbidden emission lines of all species are centrally peaked with similar kinematic distributions, indicating that the ejecta are well mixed in both SN 2024pxl and SN 2024vjm, a hallmark of pure deflagration explosion models. Radiative transfer modeling of SN 2024pxl shows good agreement with a weak deflagration of a near-Chandrasekhar-mass white dwarf, but additional IR flux is needed to match the observations, potentially attributable to a surviving remnant. Similarly, we find SN 2024vjm is also best explained by a weak deflagration model, despite the large difference in luminosity between the two supernovae. Future modeling should push to even weaker explosions and include the contribution of a bound remnant. Our observations demonstrate the diagnostic power of panchromatic spectroscopy for unveiling explosion physics in thermonuclear supernovae. 

Abstract We present SN 2023zaw—a subluminous (M_r= −16.7 mag) and rapidly evolving supernova (t_1/2,r= 4.9 days), with the lowest nickel mass (≈0.002M_⊙) measured among all stripped-envelope supernovae discovered to date. The photospheric spectra are dominated by broad Heiand Ca near-infrared emission lines with velocities of ∼10,000−12,000 km s⁻¹. The late-time spectra show prominent narrow Heiemission lines at ∼1000 km s⁻¹, indicative of interaction with He-rich circumstellar material. SN 2023zaw is located in the spiral arm of a star-forming galaxy. We perform radiation-hydrodynamical and analytical modeling of the lightcurve by fitting with a combination of shock-cooling emission and nickel decay. The progenitor has a best-fit envelope mass of ≈0.2M_☉and an envelope radius of ≈50R_⊙. The extremely low nickel mass and low ejecta mass (≈0.5M_⊙) suggest an ultrastripped SN, which originates from a mass-losing low-mass He-star (zero-age main-sequence mass < 10M_⊙) in a close binary system. This is a channel to form double neutron star systems, whose merger is detectable with LIGO. SN 2023zaw underscores the existence of a previously undiscovered population of extremely low nickel mass (<0.005M_☉) stripped-envelope supernovae, which can be explored with deep and high-cadence transient surveys. 

Abstract We present the discovery and analysis of SN 2022oqm, a Type Ic supernova (SN) detected <1 day after the explosion. The SN rises to a blue and short-lived (2 days) initial peak. Early-time spectral observations of SN 2022oqm show a hot (40,000 K) continuum with high ionization C and O absorption features at velocities of 4000 km s⁻¹, while its photospheric radius expands at 20,000 km s⁻¹, indicating a pre-existing distribution of expanding C/O material. After ∼2.5 days, both the spectrum and light curves evolve into those of a typical SN Ic, with line velocities of ∼10,000 km s⁻¹, in agreement with the evolution of the photospheric radius. The optical light curves reach a second peak att≈ 15 days. Byt= 60 days, the spectrum of SN 2022oqm becomes nearly nebular, displaying strong Caiiand [Caii] emission with no detectable [Oi], marking this event as Ca-rich. The early behavior can be explained by 10⁻³M_⊙of optically thin circumstellar material (CSM) surrounding either (1) a massive compact progenitor such as a Wolf–Rayet star, (2) a massive stripped progenitor with an extended envelope, or (3) a binary system with a white dwarf. We propose that the early-time light curve is powered by both the interaction of the ejecta with the optically thin CSM and shock cooling (in the massive star scenario). The observations can be explained by CSM that is optically thick to X-ray photons, is optically thick in the lines as seen in the spectra, and is optically thin to visible-light continuum photons that come either from downscattered X-rays or from the shock-heated ejecta. Calculations show that this scenario is self-consistent. 

Abstract We study a magnitude-limited sample of 36 broad-lined type Ic supernovae (SNe Ic-BL) from the Zwicky Transient Facility Bright Transient Survey (detected between 2018 March and 2021 August), which is the largest systematic study of SNe Ic-BL done in literature thus far. We present the light curves (LCs) for each of the SNe and analyze the shape of the LCs to derive empirical parameters, along with the explosion epochs for every event. The sample has an average absolute peak magnitude in therband of ${\overline{\mathit{M}}}_{\mathrm{r},\max }=-18.51\pm 0.15$mag. Using spectra obtained around peak light, we compute expansion velocities from the Feii5169 Å line for each event with high enough signal-to-noise ratio spectra, and find an average value of $\overline{v_{\mathrm{ph}}}=16,100\pm 1100$km s⁻¹. We also compute bolometric LCs, study the blackbody temperature and radii evolution over time, and derive the explosion properties of the SNe. The explosion properties of the sample have average values of ${\overline{\mathit{M}}}_{\mathrm{Ni}}={0.37}_{-0.06}^{+0.08}{\mathit{M}}_{\odot }$, ${\overline{\mathit{M}}}_{\mathrm{ej}}={2.45}_{-0.41}^{+0.47}{\mathit{M}}_{\odot }$, and ${\overline{\mathit{E}}}_{\mathrm{K}}=\left({4.02}_{-1.00}^{+1.37}\right)\times {10}^{51}$erg. Thirteen events have radio observations from the Very Large Array, with eight detections and five non-detections. We find that the populations that have radio detections and radio non-detections are indistinct from one another with respect to their optically inferred explosion properties, and there are no statistically significant correlations present between the events’ radio luminosities and optically inferred explosion properties. This provides evidence that the explosion properties derived from optical data alone cannot give inferences about the radio properties of SNe Ic-BL and likely their relativistic jet formation mechanisms. 

Abstract Eruptive mass loss of massive stars prior to supernova (SN) explosion is key to understanding their evolution and end fate. An observational signature of pre-SN mass loss is the detection of an early, short-lived peak prior to the radioactive-powered peak in the lightcurve of the SN. This is usually attributed to the SN shock passing through an extended envelope or circumstellar medium. Such an early peak is common for double-peaked Type IIb SNe with an extended hydrogen envelope but uncommon for normal Type Ibc SNe with very compact progenitors. In this paper, we systematically study a sample of 14 double-peaked Type Ibc SNe out of 475 Type Ibc SNe detected by the Zwicky Transient Facility. The rate of these events is ∼3%–9% of Type Ibc SNe. A strong correlation is seen between the peak brightness of the first and the second peak. We perform a holistic analysis of this sample’s photometric and spectroscopic properties. We find that six SNe have ejecta mass less than 1.5M_⊙. Based on the nebular spectra and lightcurve properties, we estimate that the progenitor masses for these are less than ∼12M_⊙. The rest have an ejecta mass >2.4M_⊙and a higher progenitor mass. This sample suggests that the SNe with low progenitor masses undergo late-time binary mass transfer. Meanwhile, the SNe with higher progenitor masses are consistent with wave-driven mass loss or pulsation-pair instability-driven mass-loss simulations. 

We present photometric and spectroscopic observations of the Type IIn supernova SN 2019zrk (also known as ZTF 20aacbyec). The SN shows a > 100 day precursor, with a slow rise, followed by a rapid rise to M  ≈ −19.2 in the r and g bands. The post-peak light-curve decline is well fit with an exponential decay with a timescale of ∼39 days, but it shows prominent undulations, with an amplitude of ∼1 mag. Both the light curve and spectra are dominated by an interaction with a dense circumstellar medium (CSM), probably from previous mass ejections. The spectra evolve from a scattering-dominated Type IIn spectrum to a spectrum with strong P-Cygni absorptions. The expansion velocity is high, ∼16 000 km s −1 , even in the last spectra. The last spectrum ∼110 days after the main eruption reveals no evidence for advanced nucleosynthesis. From analysis of the spectra and light curves, we estimate the mass-loss rate to be ∼4 × 10 −2   M ⊙ yr −1 for a CSM velocity of 100 km s −1 , and a CSM mass of 1  M ⊙ . We find strong similarities for both the precursor, general light curve, and spectral evolution with SN 2009ip and similar SNe, although SN 2019zrk displays a brighter peak magnitude. Different scenarios for the nature of the 09ip-class of SNe, based on pulsational pair instability eruptions, wave heating, and mergers, are discussed.